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SUBCRITICAL BEHAVIOR FOR QUASI-PERIODIC 
SCHRODINGER COCYCLES WITH TRIGONOMETRIC 

POTENTIALS 

C. A. MARX, L. H. SHOU, J. L. WELLENS 


Abstract. We give a criterion implying subcritical behavior for quasi-periodic 
Schrddinger operators where the potential sampling function is given by a 
trigonometric polynomial. Subcritical behavior, in the sense of Avila’s global 
theory, is known to imply purely absolutely continuous spectrum for all irra¬ 
tional frequencies and all phases. 


1. Introduction 


Let T := M/Z and fix an irrational number a, subsequently referred to as the 
frequency. Evaluating an analytic function t : T —)■ M along the trajectories 
of rotation by a with varying starting point a; G T determines a quasi-periodic 
Schrddinger operator, 


(1.1) ■= i>n-l + V’n-hl + + na)'lpn 


For every realization of the phase x 
P{Z). 


E T, (1.1) is bounded and self-adjoint on 


In physics, quasi-periodic Schrddinger operators describe the conductivity of 
electrons in a two-dimensional crystal layer subject to an external magnetic held 
of hux a acting perpendicular to the lattice plane. In this context the poten¬ 
tial sampling function v is usually a trigonometric polynomial, which, through its 
Fourier coefficients coupling constants”), carries information about the material 
properties of the crystal. The most well-known example is the almost Mathieu op¬ 
erator (AMO), in physics also known as Harper’s model, where v{x) = 2A cos(27rx) 
and A > 0. 

An interesting phenomenon encountered for quasi-periodic Schrddinger opera¬ 
tors are metal insulator transitions. Depending on the coupling constants, pres¬ 
ence of the external magnetic held may enhance or deplete the conductivity in 
the crystal. The prototype is the AMO where the spectral properties pass from 
purely absolutely continuous (ac) spectrum for A < 1 (“subcritical regime”) to 
pure point spectrum with exponentially localized eigenfunctions when A > 1 (“su¬ 
percritical regime”); the transition is marked by a “critical point” at A = 1 where 
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the spectrum is purely singular continuous (sc). For a review of the results known 
for the AMO we refer to [211 EH], 

A dynamical measure for such transitions is given by the Lyapunov exponent 
(LE), in physics interpreted as an inverse localization length, which quantihes the 
averaged asymptotics of the solutions to the time-independent Schrodinger equa¬ 
tion. Whereas positivity of the LE is heuristically associated with localization, 
zero LE is interpreted to indicate delocalization. 

Solutions to the time-independent Schrodinger equation are obtained most con¬ 
veniently in dynamical systems terms. Given an initial condition , 

n—step transfer matrices B^{x;a) allow to iteratively generate solutions of 
Hv{x);a'f’ = Exf over via 


( 1 . 2 ) 


f 

V^n-1 



0 

JJ B^{x + ja), 

j=n-l 


where 
,3) 


The Schrodinger cocycle {a,B^), a dynamical system on Tx dehned by 
{x,v) H->■ (x -|- a, B^{x)v), captures the iterative scheme in (1.2) in a compact 
way. In particular, the Lyapunov exponent of a quasi-periodic Schrodinger oper¬ 
ator is dehned by 


(1.4) 


L{a,B^):= lim — [ \og\\B^{x;a)\\dx 

n—>-oo n 


Since B^ G 517(2, M), the LE of Schrbdinger operators is always non-negative. 

In his seminal work titled Global theory of one-frequency operator^' [T], A. 
Avila introduces a framework that allows to appropriately generalize the metal- 
insulator transition observed for the AMO to arbitrary analytic potentials v. Re¬ 
lying on the analyticity of v, he considers the LE of the cocycle {a,B^{. ie)) 
obtained by complexifying the phase in (1.4); we will refer to this as the complex¬ 
ified LE and denote it by L(e; E). 

Characterized by the behavior of the complexihed LE about e = 0, Avila decom¬ 
poses the spectrum S into three mutually disjoint sets: supercritical, subcritical 
and critical energies. An energy is classihed as supercritical if the complexihed 
LE vanishes in a neighborhood of e = 0, as supercritical if the LE is positive at 
e = 0, and as critical if the LE is zero at e = 0 but not subcritical. These three 
possible situations for hxed energy E E E are shown schematically in Fig. [T] 

The supercritical regime just recovers the set of positive LE. It is the set of zero 
LE, however, for which Avila’s decomposition yields additional insight unavailable 
prior to his global theory. Whereas from Kotani-Simon theory it has been known 
that the set of zero LE forms a Lebesgue essential support for the absolutely 
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L{e;E) L{e-,E) L{e;E) 





(a) subcritical behavior (b) critical behavior (c) supercritical behavior 

Figure 1. Behavior for energies in the spectrum, corresponding to 
local behavior (at £ = 0) of the complexihed LE. 


continuous (ac) spectrum [371 EU, it still leaves unaddressed sets of Lebesgue 
measure zero where the spectrum could potentially be singular. 

Further decomposing into subcritical and critical energies, enables to explic¬ 
itly separate ac spectrum from singular spectrum. Here, as common, singular 
spectrum is dehned as the union of singular continuous and pure point spectrum. 
By the almost reducibility theorem 012] the spectrum is purely ac on the set of 
subcritical energies. It is purely singulaij^on the set of critical energies, a conse¬ 
quence of the dynamical dichotomy of Avila-Fayad-Krikorian 0. Notably, both 
these spectral results for (1.1) hold for all phases and all irrational frequencies. 

In this paper we focus on the set of zero LE and establish a sufficient criterion 
for subcritical behavior if the potential n is a real trigonometric polynomial. 


M 

(1.5) v{x) = 2 (a„ cos(27rnx) -|-sin(27ma:)) , IomI + |&m| > 0 . 

n=l 


Here, we may assume absence of a constant term which would only result in a 
shift of the spectrum. 

We mention that detecting critical energies is in principle much more delicate, 
since in contrast to both sub- and supercriticality, criticality is not stable w.r.t. 
perturbations in a and n [T]. In fact, Avila shows that small perturbations in the 
Fourier-coefficients of v destroy critical behavior which allows to prove that for 
a measure theoretically typical (=prevalent) potential v in the analytic category, 
the set of critical energies is empty [1]! 

Our sufficient criterion relies on quantifying the asymptotics of the complexihed 
LE, L(e;E), as |e| —>■ oo, building on earlier ideas (“method of almost constant 
cocycles”) which allowed to determine L(e; E) for the AMO and extended Harper’s 


^Notice that we do not claim that the spectrum on the set of critical energies is purely 
singular continuous but only that it is purely singular (see the definition of singular spectrum, 
given above) for all phases. Even though the former is a known conjecture |29j . so far this 
has not been proven yet. Excluding appearance of eigenvalues for all phases is a delicate and 
difficult problem, that has not even been established for the critical AMO. 
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model [27]. This is achieved by imposing a suitable largeness condition on 

( 1 . 6 ) 


m(e; E) := min \E — vix + ie) 

xST 


associated with the upper left entry of the matrix (1.3). We mention that such 
largeness conditions have played a role earlier in proving positivity [39l [20] as well 
continuity of the LE [TTl 0] . 

Since n is a trigonometric polynomial, m(e; E) —>■ +cx) as |e| —)■ oo. Thus letting 
{aj;bj}^^^) denote the largest e > 0 such that m{e;E) = 2, we will 

show that 
(1.7) 


L(e; E) = log \iaM — & m | + 27rM|e| , for |e| > e^- 


As we will argue, (1.7) in particular recovers the well-known lower bound for the 


LE due to Herman 
(1.8) L(a, B^) > log \iaM - &m| • 

For this reason, we will refer to ( [1.7 ) as complex Herman formula and to as 
the Herman radius. Note that as opposed to Herman’s lower bound in (1.8), the 
complex Herman formula does depend on both E as well as on all the Fourier- 
coefficients through the Herman radius. 

Exploring properties of the complexihed LE for E ^ H, will result in the fol¬ 
lowing criterion for subcriticality, which constitutes our main result: 

Theorem 1.1. Given a quasi-periodic Schrodinger operator, a irrational, and v 
as in (1.5). 

(i) E E T, is subcritical if the Herman radius, ch = satisfies 


(1.9) 


€h < 


log \iaM - &m| 
2ti{M — d) 


where 


(a) Define the uniform Herman radius e^junif = ^Jr;unif({ai; &j}^i) to be the 


(1.10) d := gcd{n € {1,..., M} : |an| + |6n| >0} . 

Define the uniform He 
largest e > 0 such that 

m(e; E) := mjn \v{9 + ie)\ = 4 -^ 2 ^ {\ian - bn\) ■ 


M 


n=l 


( 1 . 11 ) 


All energies in the spectrum are subcritical if 

log \iaM - &m| 




27r(M — d) 


Remark 1.2. 


(i) The conditions (1.9) and (1.11) implicitly assume that \iaM — 


6m| < 1 which is no restriction since the classical Herman bound (1.8) 
implies supercritical behavior if \iaM — &m| > 1- 
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ii 


We note that proving that some Eq is subcritical, implies existence of some 
ac spectrum about Eq (see also Sec. 5.2). 


We mention that both eu and en-unii can be estimated easily through polynomial 
root bounds; we discuss this in Sec. 

While the general theory is well developed, the framework of Avila’s global 
theory has only been employed explicitly to the AMO and a generalization known 
as extended Harper’s model [271 [7]. Physically interesting models include a special 
case of (1.5) with all = 0 (e.g. [231 EHl CTl HE]), also known as generalized 

focus on positivity 


Harper’s model; here the few available rigorous results 
of the LE. In light of proving subcritical behavior, we mention a related result 
on purely ac spectrum for potentials of the form v{x) = Xf{x) where / is a real 
analytic function and A G M. Here, Bourgain-Jitomirskaya prove existence of 
Ao = Xo{f) such that for all |A| < Aq, the spectrum of ( |1.1 ) is purely ac for a.e. 
a: G T if a is Diophantine na. As this result is based on proving localization 
for the dual operator, it is bound to impose arithmetic conditions on both the 
frequency and the phases. Establishing subcritical behavior, however, implies 
results on ac spectrum irrespective of such arithmetic conditions. 

We structure the paper as follows: Sec. proves the main result, Thm. 


1.1 


based on the complex Herman formula. As discussed there, the latter is an expres¬ 
sion of asymptotic uniform hyperbolicity of the Schrodinger cocycle which is quan- 
tihed by the Herman radius. The key ingredient here is Proposition |2.1| asserting 
that the Schrodinger cocycle is uniformly hyperbolic whenever m(e;E) > 2. 

Sec. [^contains a dynamical proof of Proposition |2.1 based on verifying a cone 
condition; the latter also allows to extract further estimates of the complexihed LE 
(Proposition |3.1[ ), thereby amending the result in Proposition |2^ The dynamical 
approach of Sec. [^ is contrasted with a spectral theoretical proof of Proposition 
2.1, which in particular sheds a light on the spectral theoretic meaning of the lower 
bound “2” in the largeness condition on m(e; E); as explained there, complexifying 
the phase leads to deformation of the spectrum of H^(x);a, thereby pushing a given 
energy E into the resolvent set if m(e; E) > 2. 

In Sec. we present various applications of our main result, Thm. 


1.1 to 


models of physical interest, among them to the generalized Harper model. Here, 
estimating the Herman radius is shown to be reduced to bounds on the largest 
positive root of a real polynomial, the latter of which are well explored in the 
literature. 


It is natural try to extend Theorem 1.1 to Jacobi operators, which generalize 
quasi-periodic Schrodinger operators by introducing an additional trigonometric 
polynomial c{x) whose evaluation modihes the discrete Laplacian, see (6.1). The 


extension is not immediate and is discussed in Sec. It leads to distinguishing 
three cases depending on the relative degree of c(x) and v{x). 

We conclude the paper with some remarks on how one could use the ideas from 
Sec. I^to obtain conclusions about supercritical hehavioi (i.e. positivity of the LE) 
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for quasi-periodic Schrodinger operators. Recently we learned that Jitomirskaya- 
Lin obtained the hrst qnantitive resnlts on positivity of the LE for the potential 
v{x) = 2(Ai cos(27ra:) + A 2 cos(47ra;)) which go beyond the classical Herman bonnd 
in ( 1.8[ ) [26] , Even thongh our method yields a obtain a general lower bound 
for the LE improving Herman (see Proposition |7.1[ ), it unfortunately proved to be 
difficult to extract quantitative results for a concrete potential. On the other hand, 


the lower bound in Proposition 7.1 can be easily analyzed numerically, thereby 
giving rise at least to a simple numerical scheme to test for super-criticality. 

1.1. Some notation. As common, for 5 > 0 and X = M, C,M 2 (C), the space 
C^(T, X) denotes the X-valued functions on T with holomorphic extension to a 
neighborhood of |Im; 2 | < 5, 5 > 0 equipped with supremum norm. To obtain 
statements independent of 6, we consider C‘^{T;X) := U 5 >oC 5 (T, X) with the 
inductive limit topology induced by ||.|| 5 . In this topology, convergence of a se¬ 
quence /n —)■ / is equivalent to existence of some 5 > 0 such that /„ G C^(T, X) 
eventually and \\fn — fWs —?• 0 as n —)■ cx). 


2. The complex Herman eormula 


Fix E e 


The complex Herman formula (1.7) rests on the basic observa¬ 


tion that if V is given by (1.5), the upper left corner in (1.3) will dominate the 
Schrodinger cocycle as |e| —?• 00 . Specihcally, complexifying the phase and taking 
out the dominating term yields 


( 2 . 1 ) 


B'^lx + u)= {au- — 


bu 

i 


^27rMe^—2niMx 


■1 + 0 ( 1 ) 0 ( 1 )' 

0 ( 1 ) 0 

uniformly in a; G T, as e —)■ -|-cxd. 

Thus the cocycle dynamics is asymptotically determined by the almost constant 
cocycle, ~ ))• Since L(q;,(A)^q)) = 0, continuity of the LE in the 

analytic category [27| II] implies 

(2.2) L(e; E) = log \iaM — ^m| + 27rM|e| -|- o(l) , as e —)■ cxd . 

On other hand, from Avila’s global theory the analytic properties of the com- 
plexhed LE are well understood: 

Theorem 2.1 ([I]). Let a irrational and E G M. Then, e h->- L{e]E) is convex, 
even, non-negative, and piecewise linear with right-derivatives satisfying 


(2.3) 


uj{e-, E) := —D^L{e]E) G Z 
zn 


Remark 2.2. The quantity dehned in (2.3) is known as the acceleration. Following 


the proof of (2.3) given in [T] actually shows that if v is l/d-periodic for some 
d E N, then u{e] E) G dZ. In particular for v of the form (1.5), the least positive 


value of the acceleration is d dehned in (1.10). As we will see, this accounts for 


the appearance of d in Theorem 1.1 
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From Theorem |2.1| we thus conclude existence of some Cq > 0 such that 
(2.4) L(e; E) = log \iaM — + 27 rM|e| , for all |e| > cq. 


Note that by convexity, the asymptotic formula (2.4) automatically implies a 


global lower bound, which, letting e = 0, recovers the original Herman bound in 


( 1 . 8 ). 


We mention that above argument was hrst used in m to study extended 
Harper’s model, a Jacobi operator generalizing the AMO. There, as a result of 


(2.3), the limited values of the acceleration allowed to extrapolate the asymptotics 
to obtain an expression for L{e] E) valid for all e G M. Using Remark 2.2, this 


analysis has an immediate extension to quasi-periodic Schrodinger operators if v 


in (1.5) has only one non-vanishing term, in which case for all U G S, 


(2.5) L(e; E) = maxjlog \iaM — &m| + 27rM|e|; 0} , all e G M. 

In particular, the situation is analogous to the AMO, i.e. all energies are subcriti- 
cal for \iaM — bM\ < 1, critical if |mM —&m| = 1, and supercritical if |faM —&m| > 1- 


For more general v, the simple idea underlying Theorem 1.1 is to gain addi¬ 
tional information about the complexihed LE by quantifying when the asymptotic 


formula in (2.4) holds or, put equivalently, when e h-)■ L(e; E) is eventually linear. 

To this end, we take advantage of a key result in [1] which characterizes the 
linear and positive segments of e h-L( e; E) by uniform hyperbolicity of the 
Schrodinger cocycle. The following provides a sufficient criterion for uniform 
hyperbolicity, therefore helps to identify linear pieces in the complexihed LE. We 


recall the dehnition of the auxiliary function m(e; E) in (1.6). 


Proposition 2.1. Let v G C‘^(T;M), a G T irrational, and E 
Schrodinger cocycle is uniformly hyperbolic whenever m{e; E) > 2. 


The 


Remark 2.3. The lower bound “2” of m{e-, E) is optimal in general. For instance 
if n = 0 and E = ±2, (a, B^) cannot be uniformly hyperbolic since S = [—2, 2] 
and uniform hyperbolicity is known to be an open property which, by Johnson’s 
theorem [30], cannot occur on the spectrum. The optimality will also follow 
directly from the proof given in Sec. see Remark 13 

We postpone the proof of Proposition 2.1 for now (see Sec. [^and|^ and rather 
turn to showing how it implies our main result. Theorem 1.1[ 


Proof of Theorem 1J_. First observe that for e outside the radius of zeros of 
{E — v), the minimum modulus principle implies 


( 2 . 6 ) 


m(e:E)= min \E 

\lm{z)\>e 


V Z 


whence m(e; E) increases strictly for e outside the radius of zeros of {E — v). In 
particular, the properties of e i—)■ L(e; E) stated in Proposition |2.1| imply that (2.4) 


holds with eo replaced by the Herman radius, ch = ^h{,E] {a^; 6^}“^), introduced 

in Sec. [H 
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To prove Theorem 1.1 (i), considering the contrapositive, if G S(«) is not 


snbcritical, Remark 2.2 yields oj{e = 0]E) > d. Then, the complex Herman 
formnla and convexity of L{e;E) yield the ■upper bonnd, 

0 < L(e; E) < log \iaM — &m| + ^TiMen + 27id{e — en) , 0 < e < e^ • 

In particnlar, 

log \iaM — &m| + 2TTeH{M — d) > 0, 

which is eqnivalent to 

log \iaM - bul 


> 


27r(M — d) 


Thns if en < 


logliOM - &m| 
2ti{M — d) 


, E mnst be snbcritical. 


Finally, we obtain the nniform criterion for snbcritical behavior in Theorem 1.1 
(ii), estimating the spectral radins of Tr„( 3 ;);o from above by 2 + 2 — &n|); 

this allows to eliminate the energy dependence by replacing en with the nniform 
Herman radins (ej^junif) dehned in (1.1). □ 


Theorem [lT has thns been rednced to proving Proposition 2T We shall give 
two proofs, one dynamical and one spectral theoretical, the snbjects of the Sec. 
and|^ respectively. 

3. Asymptotic domination 

We recall that factoring the Schrodinger cocycle according to 

(3,1) 

1 ‘ 


B^{x + ie) = {E — v{x + ie)) 


E—v{x-\-i€) 


) =; {E-v{x + ie))D,ix) , 


shows that its asymptotic dynamics is determined by (a, D^), which as e —>■ +oo, 
is nniformly close to the constant cocycle (a, (q o))- Trivially, the latter indnces 
the invariant splitting = ((g)) © (O')) snch that the dynamics in one invariant 
snbspace dominates the other. In this section, we prove Proposition |2. l| by showing 
that these dynamical featnres are in fact already present once m(e; A) > 2. 

To captnre these dynamical featnres precisely, we recall the following terminol¬ 
ogy from partially hyperbolic dynamics. Given D G C‘^(T;M 2 (C)) and a G T, a 
cocycle (a, D) is said to indnce a dominated splitting (also “nniform domination;” 
write {a, D) G VS) if there exists a continnous, nontrivial splitting 
and N eN satisfying 

(i) (a, T))-invariance, i.e. DN{x;a)Ex^ C for 1 < j < 2, 

(ii) for all Vj G Sx^ \ {0}, 1 < j < 2, one has 
||T)Ar(a;;Q!))ni|| ^ ||T)Af(a;;a)n2|| 


> 


+1 


+2 


(3.2) 
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Here, as in (1.2), we denote Di^{x]a) := _^D{x + ja). For obvious reasons, 

will refer to^^^ as the dominating section and to as the minoring section. 

Clearly, the condition (a, D) G VS is equivalent to some iterate of (a, D) being 
continuously conjugate to a diagonal cocycle where one diagonal entry uniformly 
dominates the other, i.e. 3N G N and C G C(T,GL(2,C)) such that 

''Ai(a;) 0 

0 X2{x), 

with Ai, A 2 G C(T, C) such that for all x G T, 

(3.4) |A,(x)| > |A2 (i)| . 


(3,3) 


C{x + Na) ^D]sf{x)C{x) = 


We mention that for analytic cocycles it is well known that analyticity is inherited 
by the invariant splitting, which in turn gives rise to analyticity of the conjugacy, 
see e.g. [H Theorem 6.1]. 

Obviously, (a, (0 0 )) ^ For a cocycle {a,D) with | detZl| = 1, VS reduces 


to the notion of uniform hyperbolicity (UT-L), in which case (3.4) simplihes to 


(3.5) 


inf |Ai(x)| > 1 


Since VS is the appropriate notion for the non-invertible cocycle (a, ( q q )), it will 
however be more convenient in this section to work with the latter. From the 
factorization in (3.1) it is clear that 


(3.6) 


a, Bf) G UH ^ (a, D^) G VS , 


whence the proof of Proposition |2.1| is reduced to showing that 
(3.7) {a, Df) G VS whenever m(e; E) > 2 . 


It is well known that VS is an open property in T x C‘^(T;M 2 (C)) [36l US], in 
particular, {a, Dfj G VS once m(e; E) is sufficiently large. The point here is to 
guantify the neighborhood of stability for VS about (a, (0 [])), which will result 
in Proposition 2.1] This will be done by verifying a cone condition^ a well known 
strategy to detect presence of a dominated splitting. 

It will be useful to work in the projective plane PC^ which we identify with 
C := C U {c)o} via (ni,n 2 ) so that D = (“ ^) G M 2 (C) acts on C \ kerU as 
the fractional linear transformation 

(3.8) z,.2:=i±^. 

a + hz 


Given a cocycle (a,/!), a coneheld for {a,D) is an open subset f/ C T x PC^ 
of the form Ua;eT{2^} x Ux such that, for all x E T, Ux is non-empty, properly 
contained in PC^, and Ux H keriA(x) = 0. A coneheld U = U3;gT{3;} x Ux for 
(a, D) is said to satisfy a cone condition if there exists G N such that for every 
X G T, one has that D]^{a]x) ■ Ux C Ux+Na- It is known (see e.g. 001) that 
verifying a cone condition implies VS. 
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Using (3.8), the proof of Proposition 2.1 is hence reduced to the following simple 
contraction Lemma: 

Lemma 3.1. For e > 0 and 0 < <5 < 1, consider the class of matrices 

{i-sy 


(3.9) 


(i^)eM 2 (C): \d\<6- \b\<e; |c|e < 


Then for each D G Se, 5 , the map Fo ■ Br(0) —)■ Br(0), Fd{z) := D ■ z is a 
contraction, where r = r{e,6) := 


Remark 3.2. The conditions on | 6 |, |c| in Lemma 3.1 are in general optimal. For 
instance, direct computation shows that if H Fd is not a contraction 


for any r > 0, cf. also Remark 2.3 


Proof. Let D G Write r^ := f, ioi some 7 G (0,1) to be determined later. If 
\z\ < r.y, then 


(3.10) 


Fn(z)l = 


|c + dz\ ^ 1 


\1 + bz\ 


7 


6'j 

c + — 
e 


Thus, Fd maps Bj. (0) to itself, if the parameters satisfy 


(3.11) 


«kl < 7((1 -~l)-S) ■ 


For fixed 6 G [0,1), the right ha nd sid e of (3.11) is maximized when 7 = 
=: 7 *, so that the condition in (3.11) becomes 

(1-5)^ 


2 

(3.12) 


e c < 


On the other hand, one has 
(3.13) 


sup 


2,loeBr.y (0) 


{Foiz) — Fd\ 

»\ 

) k 



\d — bc\ 5 + elcl 

= sup \Fd{z) \ = sup T-- :— rr < 


Z&Br.y ( 0 ) 


zeBr.y{ 0 ) 


|1 + bz \'^ (1 — 7)2 


Thus, for Fd to be a contraction on (0) it suffices to have 

(3.14) (5 + e|c| < (l- 7 )^ 
which, taking 7 = 7 *, becomes 

(3.15) e|c| < , 

in agreement with the definition of iS^^. 


□ 


Finally, we mention that Lemma 3.1 allows to extract an estimate for the com- 


plexihed Lyapunov exponent. To this end, we note that the proof of Lemma 3.1 


























11 


shows that if D = (J satishes \b\, |c| < 17 < —, then Fd is a contraction on 
i?i( 0 ) with contraction constant 

Proposition 3.1 (Proposition |2.1| amended). Letv E a G T irrational, 

and E eM.. The Schrodinger cocycle is uniformly hyperbolic whenever 
m(e; E) > 2, in which case 

(3.17) L{e;E)= \og\E — v{x + ie)\ dx + E , 

Jt 

where 

(3.18) log 

’ m(e;E)(m(e;E)- 2 ) 

Remark 3.3. It is clear from the arguments above that 
(3.20) S = L(a, Df) = 0 , for |e| > ch ■ 


_ a{e;E) \2 j_1 

m{e-,E)i ^ 


m{e-,Ep 


l + a(e;E)2 

a{e; E) = min <( 1 ; 


<2=<l°S^(l + Si 


+ 


m{f,E)^ 


m(e; E) 


Proposition 3.1 provides additional information outside the asymptotic regime, 


i.e. if m(ei; E) > 2 for some 0 < |ei| < en- In particular, for such ei, (3.18) 
implies S > — log(2), or 

(3.21) L{ei,E)> I \og\E — v{x + iei)\ dx — \og{2) . 

Jt 

We will return to this in Sec. 0 

Proof. Fix e > 0 such that m{e]E) > 2. First observe that if : T —t 
is any continuous lift of the dominating section S'i(a;;e) with ||w£(a;)|| = 1, the 
complexihed LE is given by 


(3.22) 


L{e]E)= / \og\\Bf{x)w^{x)\\ dx 


In particular, normalizing (^ the factorization in (3.1) yields (3.17) 
with 5 = L{a,Df) given by 


(3.23) 


= / log||F>,(x) 


^5i(^;^)y l|dx-- y^log(l + |5i(e;x)ndx . 

To estimate S, note that from the cone condition, S'i(x; e) is determined by the 
hxed point problem D^{x) ■ Si{x;e) = Si{x + a;e), which, by Lemma 3.1 has a 
solution in C‘^(T; i?i(0)) since 

(3.24) d = 0 , |6| = |c| = 


1 1 
\E — v{x + ie)\ ~ m{e]E) 2 


< 
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We use the following standard fact from Banach hxed point theory to obtain 
an upper bound for |S'i(x;e)|: 

Fact 3.1. Let Ti,T 2 ■. X ^ X be two contractions with contraction constants 
Ki,K 2 on a complete metric space {X,p). Denote by x* the (unique) fixed point 
ofTj, 1 < j <2. If 'j := sup 3 ,gj(^ p{Tix,T 2 x) < oo, then 

(3.25) P{xl,xl) < 7 min { ^ ^ 


Thus, comparing the solutions S'i(a:;e) and S'i(x,e = +oo) = 0, Fact 3.1 and 
(3.16) with 

(3.26) 7 = sup_ \Fdpx){z) - Fd,^^(x){z)\ < 

a;ST,2:gBi(0) 


m(e; E) — 1 ' 


implies |S'i(a;;e)| < a{e-,E) where a{e-, E) is given by (3.19). In particular, the 
upper bound on S in ( |3.18[ ) follows immediately from (3.23). 

On the other hand, using (3.24), we estimate 


(3.27) 






1 + |^i(e;a;)|2 


which implies the lower bound in (3.18) since on [0, x [0w(e;^)], {b,S) ^ 

^ 1 + 52 ^ is minimized at b rn{e-,E) ^^*3 s cl^e^Efi IT 

4. Prooe oe Proposition 12.11 - spectral theory approach 

In this section we present an alternative, spectral theoretic proof of Proposition 
2.1 Rather than verifying a cone condition as in Sec. we will use Weyl m- 
functions to obtain explicit expressions for the invariant splitting giving rise to 
uniform hyperbolicity of the Schrodinger cocycle. In particular, this argument 
will shed a light on the spectral theoretic meaning of the lower bound “2” in the 
largeness condition on m(e; E). We mention that many ideas in this section were 
inspired by our earlier work in 


We start by noting that complexifying the phase in (1.1) yields a discrete 
Schrodinger operator with complex potential, 

(4.1) := {{A+V^+i^)fi){n) = i>{n-l)+'fi{n+l)+v{x+na+ie)ilj{n) , 


in particular for e 7 ^ 0, (4.1) is a non self-adjoint operator on rfiL). Since both a 


and V are considered to be fixed, we will simplify notation by dropping the explicit 
dependence on the frequency and the potential. Henceforth, we write z := x + ze 
and use ^ and x + ze interchangeably. 

Denote by 6n the elements of the standard basis in /^(Z), and let be the 
orthogonal projection onto the subspaces, Spanjh^ : n > 0} and Span{5„ : n < 0}, 
respectively. Define the half-line operators, Hz± '■= H■ 
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For E in the resolvent sets p{Hz^±), we let 

(4.2) S-{z,E) := —m-{z,E) , 

(4.3) s+{z, E) :=—m^{z — a, E)~^ , 

where m±{z, E) := (5±i, (i^ 2 ,± ~ -S)~^<^±i) are the Weyl m—functions. The resol¬ 
vent identities show that m(. -t- ie,E) and hence s±(. -f ie,E) are continuous on 
T. The main result in this section is the following angle formula: 

Lemma 4.1 (Angle formula). Let G M. If e = Imz is such that m{e] E) > 2, 
then E G p{Hz,+) fl p{Hz-) fl p{Hz) and 

(4.4) \s.{z,E)-s+{z,E)\ = \{6o,{H,-E)-^6o)\-^>0 . 

Here, p{Hz) denotes the resolvent set of the operator H^. 

Lemma 4.1| shows that the angle between the invariant sections s±(. -1- ie,E) 
is (uniformly) bounded away from zero, in particular s±(. -f ie,E) give rise to a 
continuous, {a, i?f')- inva riant splitting of C^. 

We mention that (4.4) appeared earlier in for complex energ'ie.^as apposed 
to complex phases. Below-mentioned argument will show that the underlying fea¬ 
ture necessary in both cases is really that E G p{Hz), which of course is trivial 
for complex E and real phase. For real E, complexifying the phase leads to defor¬ 
mation of the spectrum of H^, pushing a given energy E into the resolvent set if 
m(e; E) is sufficiently large. 

Proof of the angle formula. First, we verify m{e-, E) > 2 implies E G p{Hx+ie). 
The idea is to simply write 

(4.5) -E) = A + 14+,, -E = (14+.. -E)[I+ (14+.. - E)-^A] , 

noticing tha t m(e; E) = min 3 ,gT \vix + ie) — E\ >2 guarantees existen ce o f (W — 
E)~^. (4.5) is really the operator analogue of the factorization in (3.1). The 
operator I + (Vx+ie — E)~^A is invertible if ||(l(j,+,, — i?)“^A|| < 1, which is satisfied 
since ||(l(r+,, — F^)“^|| < | and ||A|| = 2. Clearly, the same argument works for 
the half-line operators, showing that m{e] E) > 2 also implies E G p(Lfa;+i.,±)- In 
summary, all quantities in (4.4) are thus well-dehned. 


We next use some standard facts from the spectral theory of second order 
hnite difference operators, usually formulated for the self-adjoint (Jacobi) case 
(see e.g. ini)- Under the circumstances discussed here, everything is easily seen 


to carry over even though the operator (4.1) is not self-adjoint; for the reader’s 


convenience, we summarize the necessary facts including brief arguments in the 
following paragraph. 

Denote the matrix elements of the Green’s function by 

(4.6) G^{E, n, m) := (J^, - E)~^6m) ,n,meZ. 


^In fact, a continuity argument in |34j shows that for real phases, the angle formula extends 
to all E is in the resolvent of the full line operator Hx- 
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Explicit expressions for Gz{E, n, m) are available from the Jost solutions E), 
obtained by extending {Hx+ie,± — E)~^6±i G to satisfy the full line equation 

Hz'ijj = Eip. By construction, 'il!±{z,E) G does not have zeros, is P at ±cx), 
and unique up to multiplicative constants. These solutions provide the formula 


(4.7) 


G^(E,n,m) 


'p-iz, E, min{m, n})'^_|_(z, E] max{m, n}) 

W{^_iz,E),ij+{z,E)) 


verihed by direct computation. Here, W{f,g) = f{n)g{n + 1) — g{n)f{n + 1) is 
the Wronskian, which is n-independent (“conservation of the Wronskian”) if f,g 
are both solutions to Hz^|J = E'lfj. Similar computations for the half-line operators 
(eg. see [HI §1.2,2.!]) show that 


(4.8) 


m±{z,E) 


-'ip±{z,E-,±l) 

'ip±{z,E;0) 


Finally, using the same argument as in the self-adjoint case, conservation of the 
Wronskian and unicity of the Jost solutions up to multiplicative constants show 
existence of a{z, E) such that 


(4.9) 


'iIj±{z + a, E] n) = a(z, E)'ip{z, E; n + 1) . 


We note that (4.9) will l ater im ply invariance of s±(. + ie, E) under the action of 
the cocycle (a, B^) (see (4.14) below). 

To verify the angle formula, we hrst use (4.7) to express the right side of (4.4), 
thereby 


(4.10) 


KJo,(JJ.-E)-iJo)|- = 


Mz,E;l) ^p_{z,E;l) 


P+{Z:E;0) ^lJ_{z,E;0) 

On the other hand, recasting s±(; 2 , E) in terms of 'il’±{z, E), we compute taking 
advantage of (4.9), 

( 4 , 11 ) ,,(.. E) = = ppp-p = (£ - v{z)) - 


and 

(4.12) 


^/J+{z - a,E-,l) ^/J+{z,E;0) 


s_(z, E) = ————— = [E- v{z)) 


'4’+{z,E]0) ’ 


P-{z,E;0) 


i>-{z,E;0) 


Thus, combining (4.11) and (4.12) with (4.10), we conclude (4.4), as claimed. 

□ 


Lemma [4T] enters as the key ingredient in the spectral theoretic proof of Propo¬ 
sition 12.11 


Spectral theoretic proof of Proposition 2.1\ First observe that from (4.9), the sec¬ 
tions s_(. J-ie), s+(. J-ie) are naturally (a, i?f')-invariant: Under the identihcation 
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'ip±{z,E-l) 


= ^±( 2 : + a,E) . 


PC^ 3 [(i’i,i’ 2 )] ^ ^ ^ E-, one concludes: 

(4.13) 

BfM . .*(.,£) = Bf(x) = ( 

By the angle formula, s_(. + ie,E),s+{. + ie,E) thus induce a continuous, 
(a, i?®)-invariant splitting of PC^ expressed by the conjugacy, 

(4.14) a(x + (x)C.(i) = A,.,(x)-i) ' 

where Ai.e(x) and 

(4.15) 


CAx) = 




\/l — m+{z — a, E)m-{z, E) \jn-{z,E) 

are continuous in x. _ 

To conclude (a, Bf) G UH, from (3.5) it thus suffices to guarantee existence of 
N eN such that uniformly in a: G T, 


(4.16) 


AT-l 

n l^i;6(a^+i«)l > 1 

j=0 


(4.16) follows immediately if we establish L{a,Bf) > 0, in which case since 
'ip_{z,E) is P at —00 and E G p{Hz), Oseledets’ theorem determines 


n—1 


(4.17) 


n 


^log|Ai;,(x +ja)!L(a,5f) . 


j=0 


Unique ergodicity of irrational rotations shows that the limit in (4.17) is in fact 
uniform which yields (4.16). 

Finally, positivity of the complexihed LE follows from a the following well- 
known growth Lemma, which dates back to the work of Sorets-Spencer (Propo¬ 
sition 1 in [Sn]), see also [13], Chapter 3. More recent generalizations to higher 
dimensional cocycles appeared in (20], see Lemma 5.2 therein. 


Lemma 4.2 (Growth lemma). For n G N U {0}, let := 
\0‘n\ > p, some p > 2. Then for all n G N, 

1 .. 

(4.18) 


Clr). 1 


, where 


log 11 Yl Ml - log(h- 1) > 0 ■ 

j=n-l 


In particular, if = Bf{x -|- na) with m(e; E) > 2, then L{a, Bf) > 0. 


□ 
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As concluding remark we note that the proof of Lemma 4.1 can be adapted to 
show: 

Proposition 4.1. Fix a irrational and i? G M. Whenever e is such that {a, Bf) G 
lAT-L, E G p{Hz) and the angle formula (4-4) holds. 


Proof. Clearly, using the dominating (unstable) and minoring (stable) sections, 
(a, Bf') G UFL implies existence of linear independent solutions 'ip±{z] E) of = 
Eip decaying exponentially at respectively ±cx). Invariance of the sections implies 


that ip±{z]E) trivially satisfy (4.9). In summary, using these solutions in the 
formulae (4.7) shows that E G p{Hz). 

Dehne m±{z,E) from (4.8). We note that one may have m±{z,E) = oo, since 
zeros of 'ip±{z] E) are not excluded, however, the sections dehned in (4.2), 

(4.19) E) = . , 4 ,. E) = 


ij_{z,E;0) 


'ip+iz - a,E;l) ’ 


cannot both be 00 . Indeed, 'ip_{z,E;0) = fj+^z — a,E]l) = 0 would lead to 
'^+(z,i?;0) = 0 by (4.9), which in turn would imply zero Wronskian, thereby 
contradicting line ar independence of 'ip±{z]E). Thus the difference on the left 
hand side of (4.4) is well-dehned in C. Now we can run through the rest of the 
argument in the proof of Lemma 4.1 to conclude (4.4). □ 


5. Some applications 


We apply Theorem to various model situations, starting with the uniform 
criterion in part (ii) of the theorem. 

5.1. Subcriticality uniformly on the spectrum. To obtain an estimate for 
the uniform Herman radius, write (|1.5[) in complex form. 


(5.1) 


v{.x) = ^ , 

i<bl<Ar 


b' 

Aj = Qjj -|- —r , J ^ 0 , A_j = Aj . 


where ej{x) = and 
(5.2) 

Then, 

\v{x + ie)\ > \XMe-M{z) + \MeM{z)\ 


(5.3) 


M 


M-1 
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As e —)■ +CXD, the right-most side of (5.3) will eventually be positive, in particular 
letting y := eH-,nn\i can be estimated from above by the largest positive root 
Rp of the polynomial 

( M 

4 + 3^|Aj| 
i=i 

We note that p{y) has a unique positive root (Decartes’ rule of signs) an d, si nce 
p(l) < 0, necessarily Rp > 1. Thus, en-umf < Theorem 

imply: 

Proposition 5.1. All energies in the spectrum are subcritical if 
(5.5) \\m\^/^^-'^'^Rp < 1 , 


1.1 


11 


where Rp is the largest positive root of the polynomial p{y) defined in (S.f). 

Identifying subcritical behavior hence reduces to hnding Rp. 

As a hrst example, we consider the simplest nontrivial generalization of the 
AMO, letting M = 2 in (5.1). In this case, we can solve for Rp exactly, giving 

|A.I 


(5.6) 


2|A, 


+ D| \ I + I6IA2I -|- I2IA1IIA2I -|- I2IA2 

^1^2 I 


The condition in Proposition |5.1| thus yields subcritical behavior on all of the 
spectrum if 

(5.7) |Ai|-|-4|A2|-l-3|A2p-l-3|Ai||A2|<l, 

which we illustrate in Fig. 



Figure 2. The region of sub-criticality for M = 2 determined in 
(5.7). Note that the behavior captures the optimal boundary-value 
properties predicted by Fact 5.1 as A 2 —)■ 0; in this limit, |Ai| —)■ 1“. 


More generally, several articles in the physics literature consider the potential 
(e.g. |23lEHl[ni[IH]) 

(5.8) 


v(x) - 2Re{Aiei(x) + . 
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For bi = Bm = 0, this is known as generalized Harpers model, interesting also dne 
to its relation to the qnantized Hall effect in three dimensions |23l 135] . 

In this case, npper bonnds for the largest positive root of polynomials can be 
used to estimate Rp, for instance: 

Theorem 5.1 (§tefanescu [IQ]). Let p{x) = x'^ — bix'^~'^^ — ... — bkx'^~'^’^ + 
where 6 i,..., 6 *, > 0 and aj > 0 for j ^ {mi,..., nik}. Then, 

(5.9) Hi = max{(fc 6 i)i/"^b ..., 

forms an upper bound for the positive roots of p. 


Applying Theorem 5.1 to estimate Rp for the potential in (5.8), we conclude 
from Proposition 5T that all energies in the spectrum are subcritical if 

(5.10) 2|Ai| < 1 and + 6 |Ai||AMr/^^“^^ + < 1 • 


5.1.1. Limiting behavior for generalized Harper’s model. Finally, we analyze the 
limiting behavior produced by Proposition |5.1| for the potential in (5.8). Observe 
that for M = 2, (5.7) and Figure explicitly show that as IA 2 I —)■ O’*", the region 
of subcriticality approaches |Ai| —?• I"*", as expected from the spectral properties 
of the AMO. Indeed, the same behavior follows more generally from Proposition 


5.1 which we quantify in: 


Claim 5.1. For all |Ai| < 1 , there is > 0 so that for all 0 < \Xm\ < L some 
n = fi:(|AM|), Proposition 5.1 guarantees subcritical behavior on all of the spectrum 
whenever IA 2 I,..., |Am-i| < Specifically, writing 0 < |Ai| = 1 — 61 — S 2 , for 
61,62 > 0 , one can take 

/ r ^ ]\4 —1 

(5.11) /i = ' ^ 


and for M >2, K= 


PmI 

2(M-2) 


4 + 3M 

M-2 


Proof. Write |Ai| = 1 — (5i — ^2 for 61,62 > 0. Note that from (5.11) both p,, k < 1. 
Letting 


(5.12) 


A := (M — 2) max |A,| < k(M — 2) , 

2<|i|<M-2 


M-U 


we estimate 

(5.13) p{y) > IXmIv^ -y + { 6 iy - (4 + 3M)) + {62y - Xy^' 
Lemma 5.2. Let yo := ■ Then y > yo implies p{y) > 0. 


Proof. Let y > yo. First note that (5.11) implies y > so that 

(5.14) { 6 iy - {4 + 3M)) > 0 . 
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Denote by zq the unique positive root of |Am||/^ —(^ —2)i/^ ^ — For M = 2, 
^0 = 2 / 0 - Using Theorem 5T, for M > 2 we estimate 

2(M - 2) 


(5.15) 


^0 < 


|A 


M\ 


We distinguish the cases y < Zq and y > Zq. y > Zq, then \Xm]^ — Xy^~^ — 
y > IXnly^ — {M — 2)y^~^ — i/ > 0, so that the right side of (5.13) is positive, 
whence the claim in the lemma is satisfied. 

\i y < Zq, then for M > 2, the definition of k and A as well as (5.15) imply, 

M-2'' 


(5.16) 


{ 62 y - \y^ ^) > ^22/ I 1 


/ zo I Am I 
V2(M-2) 


Thus for y>yo= , (5.13) yields 


Piy) > |Am||/^ -y+ {Siy - (4 + 3M)) + {62y - Xy^^^ > {XmW^^ - 2 / > 0 , 


M-l\ 


> 0 . 




>0 


— 

>0 


as claimed. 


□ 


We conclude from Lemma 5. ^that Rp < I/O, whence 


(5.17) 


Rp\X 


M 


1 /(M- 1 ) 


< 1 , 


which in particular implies the condition (5.5) in Claim 5.1 


□ 


We note that qualitatively this behavior is expected from the point of view of 
Avila’s global theory; it is an immediate consequence of upper-semicontinuity of 
the acceleration in the cocycle [1]: 


Fact 5.1. Given a quasi-periodic Schrodinger operator Hy(^xy,ao with analytic po¬ 
tential vq and irrational frequency ao? suppose that all energies in the spectrum of 
Hvoix)-,ao o,re suhcritical. Then the same is true for all {v,a) G C^(T;M) x M \ Q 
in some open neighborhood of (no,ao). 

Proof. By Theorem m the acceleration for e > 0 can only attain non-negative 
integer values if a is irrational. Hence, upper semi-continuity of the acceleration 
and compactness of the spectrum ensure existence of some open neighborhood 
of (ao,no) in C‘^(T;M) x M \ Q, such that u(e = 0,E) = 0 for all E in the 
spectrum. Here, we also use that the spectrum of quasi-periodic Schrodinger 
operators depends continuously on {a,v) in the Hausdorff metric [S]. □ 


We note that in spectral theoretic terms, Fact |5.1| states that purely ac spectrum 
is stable w.r.t. perturbations in {a,v). 

Unfortunately, as Ai —)■ 0+, Proposition ^ does not capture the optimal be¬ 
havior predicted by Fact 5.1 Indeed, since fh{e;E) depends continuously on the 
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Aj’s, we set Ai = 0 in (5.8) and explicitly solve for ej^junif to find 

IjM 


(5.18) = 

which, if eH;unif < 

(5.19) 



+ 1 + 


|A 


+ 


M\ 


|A 


+ 2 


> 


M 


|A 


M 


1/M 


27r(M-l) 


imposes the restriction 
1 


|A 


M 


< 


A 

4^ 


1 

< - 
- 4 


5.2. Energy dependence. To illustrate the application of Theorem Ow , we 
consider odd potentials where all Oj = 0 in (1.5). In this case, it is known that 
E = 0 G S(a), indeed: 

Fact 5.2. Given a quasi-periodic Schrodinger operator such that for some Xq G T, 
v{. + Xq) is odd. Then, for every irrational a, 0 G S(q;). 

be found in [10]; for the reader’s convenience we 


can e.g. 


A proof of Fact 5.2 
give a slightly shorter, alternative argument in Appendix 

As an example, we consider the simplest case where v{x) = hi sin(27ra;) + 
^2sin(47ra;). Estimating like in (5.3), allows to bound the Herman radius for 
E = 0 from above by the largest positive root Rq of the polynomial, 

(5.20) q{y) := \b 2 \y‘^ - \bi\y - (|6i| + + 2) , 

where, as before, y = Thus, the condition (1.9) in Theorem 1.1 
that E = 0 is subcritical if < 1, which, computing Rq, yields 


fi) asserts 


(5.21) 


|&2p + \bi\\b2\ + 2I62I + \bi\ < 1 


As mentioned earlier (see Remark 1.2| (ii)), proving that E = 0 is subcritical 
implies existence of some ac spectrum centered around E = 0. This follows again 
using upper-semicontinuity of the acceleration, which since E = 0 is subcritical, 
implies that u{e = 0, E) = 0 for all E in some interval (^1,^2) containing 0. 
Hence, by the almost reducibility theorem, all spectral measures are purely ac on 
S(q;) n (El, E2). 

Fig. I depicts the region determined by ( |5.21[ ) where the operator has cor¬ 
respondingly some ac spectrum; the same hgure compares this with the region 
determined by (5.7) where all of the spectrum is purely ac continuous. 

More generally, using the same ideas for 

(5.22) v{x) = 2{bN sm{27rNx) + bM sin{2nMx)) , 1 < N < M, bM ^ 0 , 


we conclude, employing the root bound in Theorem |5.1| , that: 

Example Theorem 5.3. Consider a quasi-periodic Schrodinger operator with 


potential given by (|5.22|) 
(5.23) \b 


M 


1/M-N 


E = 0 is subcritical if 

1/(M-Af) 


max 


\b 


N\ 


f 4 -I- 2|&jv 


\ l/AT 

+ 2 


< 1 . 
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Figure 3. For the Schrodinger operator with potential v{x) = 
6i sin(27ra;) + 62 sin(47rx), the upper curve represents the boundary 


of the region determined by (5.21); below the upper curve, E = Q 


is subcriticah This is compared with the region according to (5.7) 


where all energies in the spectrum are subcritical; the lower curve 
represents the boundary of the latter. 


In particular, the operator has some ac spectrum if (5.23) is met. 


6 . Jacobi operators 


6.1. Jacobi cocycles. It is natural to try to extend Theorem 1^ to Jacobi op¬ 
erators. 


(6.1) := c{x + {u - l)a)V'n-i + c(x na)'i/^ri+l + v(x na)i’r, 


Here, the discrete Laplacian in (1.1) is modihed by evaluating the (complex) 
trigonometric polynomial 

N2 


( 6 . 2 ) 


C(X = 


k=Ni 




Jlirikx 


Ni < N2 i \^Nl^^N2 \ > 0 , 


along the trajectory x ^ x + a. As before, v is assumed to be (real) trigonometric 
polynomial of the form given in (1.5). 


A prominent example from physics is extended Harper’s model (EHM) where 
both c, n are both trigonometric polynomials of degree 1. Proposed by D. J. 
Thouless in context with the integer quantum Hall effect |12], EHM generalizes 
the AMO, allowing for a wider range of lattice geometries by permitting the 
electrons to hop to both nearest and next nearest neighboring lattice sites. Its 
spectral theory has recently been solved [7], relying on an extension of parts of 
Avila’s global theory to analytic Jacobi operators [271128]. In fact, the “method of 
almost constant cocycles” underlying the complex Herman formula was originally 
developed in [27| to hnd the complexihed LE for EHM. 
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As has been detailed in [2^ 155] . the complexihed LE for Jacobi operators is 
defined bjj^ 

(6.3) L{e; E) := L{a, Af) — I{c) , I{c) -.= f log |c(a;)|da: , 

Jt 

where L{a, Af) is the LE of the phase-complexified Jacobi cocycle induced by 


(6.4) 


Afix) : = 


E — v{x + ie) —c{x — ie — a) 
c(x + ie) 0 


Since /(c) is a constant independent of e, determining the complexified LE (6.4) 
reduces to finding the LE associated with (6.4). 

As before, letting e = 0 in ( |6.3 ) yields what is usually known as the LE of a 
Jacobi operator (6.1). We mention that alternative choices for Jacobi cocycles 
exist (see e.g. P4l 129] for details), however for what is to come, (6.4) turns out to 


be the most advantageous. 

One feature not present for Schrodinger cocycles is that {a, Af) is in general 
non-invertible, indeed, det Af(x) = c(x — ie — a)c{x + ie) which may vanish due 
to zeros of c. A Jacobi operator is called singular if c{x) has zeros on T and 
non-singular otherwise. Singularities of the cocycle often lead to interesting phe¬ 
nomena when trying to generalize results for Schrodinger operators to the Jacobi 
case, which has been explored in several recent articles |29lES SSI laiia EH 12711111 
m [23] . From a dynamical point of view, presence of singularities is accounted for 
by replacing uniform hyperbolicity {WH) with uniform domination {T>S) (recall 
Sec. [^for a dehnition) Sj- 

It was proven in [271128] (see also [29]) that Theorem |2 .1 j essentially carries over 
to Jacobi operators: 


Th eorem 6.1. The analytic properties of e L{e]E) stated in Theorem 
2.1 hold essentially unchanged with the only alteration that u{e; E) G 
For non-singular Jacobi operators one still has a;(0; E) G Z, in particular 
the smallest non-zero value of uj{0] E) is 1. 


In particular, for both singular and non-singular Jacobi operators. Fig. [T] rep¬ 
resents the three possible situations for the graph of L(e; E) in a neighborhood 
of e = 0 if E G E. From a spectral theoretic point of view the implications 
for each of these three cases for non-singular operators are the same as in the 
Schrodinger case [3129]; one thus partitions E into subcritical, supercritical, and 
critical energies. 

For singular operators, partitioning the set Z := {E : L(0; E) = 0} into 

subcritical and critical does not yield any further information, as the presence of 
zeros of c{x) a priori excludes any absolutely continuous spectrum for all a; G T 


^Replacing the complex conjugate of c by its reflection along the real line as done in the upper 
right corner of (6.41 makes (6.41 an analytic, matrix-valued function. 
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[T^ . For singular Jacobi operators, the criteria derived below will thus simply 
imply zero LE. 


6.2. Asymptotic Analysis. We now turn to generalizing Theorem 1.1 to Jacobi 
operators, and will encounter two complications: 

First, observe that the method of almost constant cocycles does not generalize 
immediately to t he J acobi case - if we naively factor out the leading terms in 
Af in analogy to (2.1), the remainder for the case deg(c) > deg(n) will in general 


approach a constant matrix with zero spectral radius (and thus LE = —cxd), which 
would yield an undetermined expression for L{a,Af) in the limit e —>■ oo. 

This problem is remedied by hrst conjugating Af by 


(6.5) 


C = I n K 


for some appropriate K E where for convenience we write w := 

Here, we call two cocycles {a, A) and {a, A') with A, A' G C‘^(M/2Z, M 2 (C)) con¬ 
jugate if A'{x) = C{x + a)~^A{x)C{x) for some C G C'^(M/2Z, GL(2, C)). Con- 
jugacies clearly preserve the LE. 

Applying the conjugacy in ( 6.5| ) will lead to consideration of cases, depending 
on the sign of 2M — {N 2 — Ni). The latter expresses the dependence of the 
asymptotics of L{e-, E) on the relative degree of c and v. We note that from a 
spectral theoretic point of view, the conjugacy in ( 6.5| ) is equivalent to a known 
unitary whose action transforms the original Jacobi operator to one where c{x) is 


replaced by c{x) 


^—2'KiKx 


see e.g. [H], Lemma (1.57) and Lemma 1.6, therein. 


The second complication is of fundamental nature. The key in Sec. was to 
quantify when the asymptotics expressed through the complex Herman formula 
holds. This was possible since Schrodinger cocycles are asymptotically (as e —)■ 00 ) 
close to the constant Doo = ( 00 ) with (a, Hoo) £ ES; quantifying the asymptotics 
then amounted to Ending the radius of stability of VS about (a, Doo)- 

Depending on the sign of 2M — {N 2 — Ni), this will not be possible for Jacobi 
operators. Indeed analyzing the case 2M — {N 2 — W) < 0 will show that the 
limiting constant cocycle lies on the boundary of VS. The asymptotic expression 
however still leads a Herman bound which, due to above remarks, is not entirely 
trivial and has not appeared before. 

We turn to the analysis of the cases, starting with 2M > {N 2 — W), where a 
criterion for sub-criticality can be obtained. Here, the analogue of (1.6) is defined 
by 


( 6 . 6 ) 


, \E — vixie)\^ 

mie; E) := mm ----r 

xeT \c[x -1- ie)c[x — le — a) 
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respectively, uniformly over E E T,, 

(^|u(a; + ie)| - 2 \f^j\ + T.f=i l^l)) 


(6.7) m(e: E) -.= min 

xGT 


\c{x + ie)c{x — ie — a) 


where Xj is given in (5.2), as earlier. 

Define the corresponding Herman radii, en and e^junif, as the largest e > 0 such 
that, respectively, m(e; E) = 4 and m(e; E) = 4. 

Theorem 6.2. Suppose 2M > {N 2 — Ni). 

(i) Then for every E eM., 

(6.8) L{a, A^) = 2'KK\e\ + log |Am| ; M |e| > ^h- 
In particular one has the Herman bound, 

(6.9) L(0;E) > loglAjul-/(c) . 

(a) For non-singular Jacobi operators, E E S(a) is subcritical whenever 

^ -loglAjul+/(c) 

^ 2MM - 1) ■ 

All of the spectrum is subcritical whenever 

^ -loglA^I+/(c) 

(6.11) eH;unif< 27r(M-l) ' 

For singular Jacobi operators, ( 6.10^ and (6.11) simply imply F(0;F) = 0. 

Proof. Let K E to be determined later. Conjugating the Jacobi cocycle by C 
in (6.5), we obtain 

( 6 . 12 ) A'^{x) := 


E — v{x -1- ie) 
Z-K^2niKac(^X -\- ic) 


-w^c{x — ie — a] 

0 


Note that taking the limit e -E +cxd is equivalent to |w| —)■ 00 . Thus, expressing 
the two off-diagonal terms of (6.12) in terms of w, 




N2 

(6.13) 

w^c{x — ie — a) 




j=Ni 

(6.14) 

yj-K^ 2 niKa^^^ -f- ie) 

N2 


j=Ni 


we see that the upper left corner in (6.12) will dominate as e —>■ J-cx), if 


N 2 + K <M , 

(6.15) iVi + F:>-M, 

which is possible since iV 2 — iVi < 2M. 
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Ari+jV2 


Now we can factor out the 


In agreement with (6.15), we pick K = ^ 

dominating term in A'^ and take e —)■ +oo, giving 

(6.16) L(a, Af) = L(a, A',) = 2i,Mi + log |A„| + L(a, (|^ 


where 

(6.17) /„(t.i) = 

(6.18) 


E — v{x + ie) 

^^^27vMe—27viMx 


1 , \fb{e,x)\ = 

l/c(e,x)| = 


cur 


ze 


a 


\XM\e^AM-K)e 

\c{x + ze)| 


|AM|e2-(^+^)^ ■ 

Since N 2 — Ni < 2M, fa, fb, fc = o(l) uniformly in x as e —)■ +cxd, hence 


(6.19) 


{a,Df) := (a, 


l + /a(e,x) /b(e,x) 

fcie,x) 0 


—)■ [a, 


1 0 
0 0 


G VS . 


Thus, from the stability statement in Lemma 3.1, we conclude {a,Df) G VS if 
e > en- 

On the other hand it is known that BIB L(e; E) is li near a nd posi tive if and 
only if {a,Af) G VS, which, combining Theorem |6 .1 1 and 6.16, yields (6.8). Note 
that as in the Schrodinger case, we use that e 1 —)■ m(e; E) increases strictly for 
e > Finally, convexity of the complexified LE implies the Herman bound, 

(0. 

To prove part (ii), we first assume that the Jacobi operator is non-singular and 
follow the same contrapositive argument as in the proof of Theorem 1^ If E G S 
is not subcritical, Ci;(e = 0; E) > 1, which by convexity would imply the upper 
bound 


( 6 . 20 ) 


0 < L(e; E) < log |Am| + 2756^ + 27r(e - eu) , 0 < e < e^- 


Notice, that by Theorem 6.1| a;(0; E) > 1 for non-singular operators if E G S is 
not subcritica0 

In particular, letting e = 0, we obtain (6.10) upon taking the contrapositive. Fi¬ 
nally, the uniform condition in (6.11) follows immediately estimating the spectral 
radius of Hc(x),v(x)-,a from above by 

( 6 . 21 ) 


2||c||oo + ItIIoo < 2 


M 


i=i 


|A, 


N 2 

.j=Ni 

Finally, if the operator is singular, we can use density of non-singular Jacobi 
operators in operator-norm topology (see Lemma 6.3, below) to extend the upper 
bound in (6.20) to the singular case, which then yields the claim in (ii): To see 
this, by the proof of Lemma 6.3, there exists a sequence (c„) of trigonometric 


'^Indeed, here we only use the “only if” direction which is essentially trivial. 

^For singular operators, the least positive value would be | which would immediately imply 
a weaker form of (6.20). The limiting argument below however allows to improve on that. 
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polynomials such that c„ —>■ c in C‘^(T), and for all n G N, c„ has upper and lower 
degrees N 2 and iVi, respectively, and has no zeros on T. Note that this in particular 
implies that the condition 2M > {N 2 — Ni) holds along the resulting sequence of 
approximating non-singular Jacobi operators, which allows application of above 
argument for the non-singular case. 

Denote the spectrum of Hc„(x),v(x)-,a by and the spectrum of Hc^^x),v{xy,a by S. 
Clearly, for all x e T, Hc^(^x),v(xy,a -t Hc^x),vixya in norm-topology, in particular, 
S in the Hausdorff metric. 

Suppose E E T, does not satisfy L{e = 0;E) = 0. Then, taking En E 
such that Eji -E E, continuity of the LE of analytic cocycles [27], implies that 
L{e = 0] En) > 0, eventually; in particular, En is not subcritical whence (6.20) 
holds eventually. Taking the limit implies that the upper bound in (6.20) holds 
for E. □ 


Lemma 6.3. The set of functions in C‘^(T) which are hounded away from zero on 
T is open and dense in C‘^(T). In particular, non-singularity of analytic Jacobi 
operators is Baire-generic in operator norm. 


Proof. Openness is clear. To show density, given 0 ^ f E C‘^(T) that has zeros on 
T, factorize f{x) = t{x) ■ q{x), where q is zero-free on T and t{x) is a trigonometric 
polynomial containing all the zeros of / on T. Then let (e„) be a real sequence 
with Cn = 0 ( 1 ), letting 


( 6 . 22 ) 


fn{x) := t{x + ien)q{x) , 


fn has no zeros on T and /„—;■/ in C‘^(T). 

We note that if / is a trigonometric polynomial, then /„ is a trigonometric 
polynomial of the same degree than /; the latter is relevant for the proof of 
Theorem |6.2| (ii). □ 


Note that in order to apply the stability Lemma |3.1| , the previous argument 
relied on m(e; E) -E J-oo as e —)■ J-oo, in particular m(e; E) = A eventually. If 
2M = {N 2 — Ni), this is not the case anymore, in fact 


(6.23) 


Lemma 


m(e; E) —)■ ^ (r -goo . 


|A 


M 


3.1 


is however still applicable if ^ < |, in which case both ch and 


eH;unif are well-defined. Hence, we conclude: 


Theorem 6.4. Suppose 2M = N 2 — Ni. 

(i) For all E eM. and some 0 < cq, 
(6.24) 

L{a, Af) = 27rM|e| -|- log ( max 




all |e| > Cq. 
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(6.25) 


(6.26) 


In particular, one has the Herman bound 

Am ± 


L(0; E) > log max 


l ± 


-/(c) 


If 


l/'-ATi/'-jVa 

|AmP 


1 

^ 4 


then ( 6 . 24 ) holds for all |e| > ch- 
(ii) Suppose (6.26) holds. 

If the operator is non-singular, then E G S(a) is subcritical whenever 


log max± 


(6.27) 


< 


Am ± -x/Aif ~ 4/iAr2/iAr^e2’^*^" 


+ I{c) 


27r(M - 1) 


All of the spectrum is subcritical whenever (6.27) holds with cr replaced 
by Ci7;unif’ 

For singular Jacobi operators, (6.26) and (6.27) imply F(0; F) = 0. 


Remark 6.5. We note the frequency dependence in the lower bound (6.25). Indeed, 
it follows from our earlier work on the LE of extended Harper’s model in 123, that 
for W = — IjW = 1 and M = 1, the frequency dependence of the asymptotics 
(6.24) persists as e —)• 0, resulting in a frequency dependence of the LE of the 

Ai ± a/A^ - 4FI/i_ie2™ 


Jacobi operator, 


(6.28) F(0;F) = max log ( max 


-J(c);0 


This is interesting, since for quasi-periodic Schrodinger operators there is no 
known example of a LE with explicit dependence on a. 


Proof. The argument follows the same steps as in the proof of Theorem 6.2 We 
again conjugate by C* = (J Jk), this time with K = M — N 2 = —Ni — M. 
The leading terms in the off-diagonal entries of (6.12) are n ow and 

^ 2 -KiKa^^^^M^ Thus, we can pull out and use Theorem 6.1 to see that for |e| 
sufficiently large, 

(6.29)F(a,Hf) = 27rM|e| + logp 2 M(M-iV 2 )a 


Am Rn 2^ 


2TviN20^ 

0 


(6.30) 


= 2TiK\e\ + log I max 


e '- Rni 

Am ± V^M~ 


which by convexity implies the lower bound in (6.25). In (6.29), we use p{M) to 
denote the spectral radius of a matrix M. 
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Moreover, if (6.26) holds, the limiting constant cocycle in (6.29) induces a 
dominated splitting, whence we conclude from Lemma 3.1 that (6.30) holds for 
|e| > en- Here, we also use that concavity of e h-)■ logm(e;£') (Hadamard’s 
three-circle theorem) and the maximum modulus principle imply that m{e]E) 

^ strictly as < e —t -|-oo. 

Part (ii) of the theorem is now obtained using identical arguments as in the 
proof of Theorem |6.2[ □ 


Lastly, we turn to the case when 2M — {N 2 — Ni) < 0. Then, the conjugacy 
mediated by C* in (6.5) does not resolve the problem of an undetermined expression 
when considering the asymptotics. 

Instead, we consider the second iterate of the Jacobi-cocycle, (2a,Af.2), where 
= Af{x -f a)Af(x) with 


Afl^ix + ie) = 


{E 


v{x)){E — v{x -|- a)) — c{x — ie)c{x) —{E — v{x + a))c{x — ie — 
c{x + a){E — v{x)) —c{x — ie — a)c{x + a) 

and L{2a, Af. 2 ) = 2L{a, Af). Then, since 2M < N 2 — Ni, we can write 


(6.31) Af.^2i^ + te) = e^AN.-N,)e^-2.^iN,-N,). l^o(l) - 


0(1) 
C’(l) “ 


as e —)■ -|-cx). 

By Theorem |6.1 we thus conclude: 

Theorem 6.6. If N 2 — Ni > 2M, then for all E G 
that 


i, there exists 0 < cq such 


(6.32) 


L{a,Af) = 'k{N 2 - A^i)|e| -f -(log I/iatJ log I/xtvJ) all |e| > cq. 


In particular, one has the Herman bound, 


(6.33) 


L{f)-,E) > ^(log|/xjVil Tloglhval) -I{c) . 

7. Some remarks on supercritical behavior 


In this hnal section, we comment briefly on how one could use ideas from Sec. 
to obtain conclusions about supercritical behavior (i.e. positivity of the LE) 
for quasi-periodic Schrodinger operators. The following relies on the estimates 
of the complexihed LE in Proposition 3.1, which, assuming existence of some Ci 
satisfying 

(7.1) 0 < ei < 6// with m{e,E) > 2 . 

will allow to extract a lower bound for (see (7.4), below), thereby im¬ 

proving on the classical Herman bound ( |1.8 ). 

Testing for existence of such ci requires estimates of the function m(e; E) outside 
the asymptotic regime which, unfortunately, we have found difficult to extract. 
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On the other hand, it is easy to solve for m(e; E) numerically, which at least gives 
rise to a simple numerical scheme to test for supercritical behavior. Below, we 
will demonstrate this for generalized Harper’s model (5.8) with M = 2. 

Assuming existence of ei satisfying (7.1), we first establish above mentioned 
improvement on the Herman bound: 

Proposition 7.1. Consider a quasi-periodic Schrodinger operator with trigono¬ 
metric potential, v{x) = = ^-j a irrational. Given 

E eM., suppose there is Ci = ei{E; satisfying (7.1). Then 

ch -e 


(7.2) 

where 


L(e; E) > log |Am| + 27rMe + 7 


(-H — Cl 


0 < e < ei 


(7.3) 7 := / log \E - v{x + iei)|dx - log( 2 ) - log |Am| - 27rMei . 


In particular, letting e = 0, one has 


(7.4) 


L{a,B^) > log|AM| + 7 - 


en 


— ei 


Remark 7.1. Using Jensen’s formula, the integral fj log \E — v{x + zei)|da; can be 
evaluated based on the zeros of the polynomial 

M 

(7.5) fE,e,{z) := Ez^ - . 

i=i 

Letting oi,..., be the zeros of fE,ei{z) (counted with multiplicity) in the com¬ 
plex unit disk D, 7 is given by 

n 

7 = ^ log 


(7.6) 


k=l 


\ak\ 


log 2 . 


Proof. Consider the line segment connecting the points (ei, L{ei;E)) and (eo, L{eo;E)). 
By convexity of e h->■ L(e; U), one necessarily has 


, ^ L{eH;E)-L{euE) 


(7.7) 
whence 

(7.8) L(e; E) > L{ei]E) 27ra(e — ei) , 0 < e < ei . 

On the other hand, one has L^ch', E) = log |Am| +27rMeH, and the lower bound 
(3.21) from Proposition |3.1| implies 

(7.9) L{ei]E)> / log lU - n(a:- 1 -iei)|dx - log(2) = 7 -Flog |Am| 27rMei , 


where we have made use of the dehnition of 7 in (7.3). 
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Thus, we can estimate 


(7.10) 


a = 


L{eH;E)-L{er,E) 
27i{eH - ei) 


< M 


7 


2TT{eH 


which, combined with (7.8), yields (7.2). 


□ 


7.1. Example (numerics). Consider n(a;) = 18 cos(27rx) + 1.6 cos(47rx), so Ai = 
9 and A 2 = 0.8. Using Mathematica, we hrst computed numerically m(e; E) = 
minj-gT \v{x + ie)|, the results of which are shown in Fig. 



Figure 4. Plot of m(e; E) for Ai = 9, A 2 = 0.8. The horizontal 
line is drawn at 4 + 2Ai + 2 A 2 = 23.6. 


At e = 0.2, minajgT \v{x + i0.2)\ ^ 24.242 > 23.6 = 4 + 2Ai + 2 A 2 , so that 
m(e; E) > 2 for all E G [—21.6, 21.6] 3 S; thus for any E in this interval, we can 
apply Proposition |7.1| with ei = 0.2. 

For a specihc example, take E = —2. The roots of f{z) = —2z^ — — 

9e2^(0'2)^-1.6e-4U0.2)^4_^_gg47r(0.2) -39.055, -0.3161, and -0.07839±3.51271i 

Only one root, z = —0.3161 is in i?i(0), so from (7.6), 

(7.11) 7 = - log(0.3161) - log 2 = 0.458. 


To apply (7.4), we hrst solve numerically solve for the Herman radius, which 


for E = —2, yields en = 0.3864. Then, the lower bound in (7.4) implies 
(7.12) L{a, B^) = L(0; E =-2) > log(0.8) + 0.458 


.3864\ 


.1864 / 


= 0.727 > 0 


In comparison, the classical Herman bound gives L(0; E) = L(q;, B^) > log |Am| ~ 
-0.223, for all E G M. 

Using Mathematica, we also sampled energies E G [—21.6,21.6] ^ S, using 
a step size of 0.001. We simplihed the computation, by computing the uniform 
Herman radius instead of en for each energy. From the proof of Proposition 


7.1 


it is clear that (7.2) also holds for en replaced by en-^unif-, since then still 
L{eH ur,\f \ E) = log I Am I + 27 re//;unifAf which was used in (7.9). The computation 
leading to Fig. [fallows to extract the uniform Herman radius, ej^^unif = 0.4142. 
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Applying the bound in (7.4), numerical computation of 7 using (7.6), results 
(7.13) 


L(0; E) > log IA 2 I + 7- 


eH;unif 


We show a plot of L~{E) in Fig. 


=: L-{E) > 0 , for all E E [-21.6,21.6] C S. 


Lower bound for L(0;E) 



Figure 5. Plot of the lower bound L (E) for the LE extracted 
from Proposition |7.1| for Ai = 9, A 2 = 0.8, and E E [—21.6,21.6] 3 

E. 


Appendix A. Proof of Fact 15.21 

We approximate a by some rational in which case the resulting discrete 
Schrodinger operator becomes g-periodic. From the theory of periodic Schrodinger 
operators, it is know that for each x G T, the spectrum a{^,x) is determined by 

the discriminant E) = tr x)| via 

(A.l) ff(^,i)=A(?i;.)-'([-2,2|). 

We claim: 

Lemma A.l. Let Xq as in Fact \5.S\ then for all E eM., 

(A.2) A(^, xo; E) = (-1)''A(^, xo; -E) . 

q q 

In particular, if q is odd, then E = 0 E cr(|,Xo). 

Proof. For simplicity denote Vj = v{xq + j^),0<j<g — 1, and, for /i G M, set 
:= Note that for every qi, 

(A.3) a, = (-1)A7 , 
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where superscript t is the matrix transpose. Thus, using the anti-symmetry of 
v{xq + .) and invariance of the trace under matrix transposition and cyclic per¬ 
mutation, we obtain 

A(?, B) = tr { (nf.-y-'" A^-v) Ay:} 

= (-l)»tr I Aty-V^ (nj‘.(cI-l)/2 A'_yyy^ A'_y} 

= (-l)«tr {a.,, (n'rd'" A_yyy) (n‘.„-„/2 A.y.v) ) 

(A.4) = (-l)'>A(E.i„i-B) . 

□ 


Finally, let (|^) be any sequence of rationals with ^ a and odd for all 


q-n 


n G M. Since the map M 3 /3 i— )■ S+(/3) := is known to be continuous 


in the Hausdorff metric at every irrational a [8], Lemma A.l implies Fact 5.2 
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